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ABSTRACT 

We investigate the dust associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC) as observed with the Herschel 
Space Observatory. N49 is unusually bright because of an interaction with a molecular cloud along its eastern edge. We have used PACS and SPIRE 
to measure the far IR flux densities of the entire SNR and of a bright region on the eastern edge of the SNR where the SNR shock is encountering 
the molecular cloud. Using these fluxes supplemented with archival data at shorter wavelengths, we estimate the dust mass associated with N49 
to be about 10 Mq. The bulk of the dust in our simple two-component model has a temperature of 20-30 K, similar to that of nearby molecular 
clouds. Unfortunately, as a result of the limited angular resolution oi Herschel at the wavelengths sampled with SPIRE, the uncertainties are fairly 
large. Assuming this estimate of the dust mass associated with the SNR is approximately correct, it is probable that most of the dust in the SNR 
arises from regions where the shock speed is too low to produce significant X-ray emission. The total amount of warm 50-60 K dust is ~0.1 or 0.4 
Mq, depending on whether the dust is modeled in terms of carbonaceous or silicate grains. This provides a firm lower limit to the amount of shock 
heated dust in N49. 

Key words. ISM: supernova remnants: SNR N49; ISM: dust; Galaxies:Magellanic Clouds; Submillimeter: ISM 



> 

oo 

in 

o 
o 



X 



1. Introduction 

N49 is a bright X-ray and optical supernova remnant (SNR) in 
the Large Magellanic Cloud (Long et al. 1981), associated with 
a soft gamma ray repeater (Cline et al 1982) and unresolved X- 
ray point source (Rothschild et al. 1994). It has an asymmet- 
ric surface-brightness distribution at essentially all wavelengths, 
which is due to an interaction with a molecular cloud that is lo- 
cated on the southeast limb of the SNR (Vancura et al. 1992; 
Banas et al 1997). Based on a simple Sedov model for the SNR, 
Hughes et al. (1998) concluded that the swept-up mass was the 
order of 200 Mq. Park et al. (2003) recently conducted a detailed 
study of its X-ray properties using Chandra, confirming that the 
SNR is relatively young, on the order of 6 000 yrs. They find that 
although most of the emission is dominated by interstellar mate- 
rial (ISM), the Chandra observations appear to show "explosion 
fragments" on the SW side of the SNR in X-rays. According to 
Bilikova et al. (2007), the X-ray temperature (-1.2x10^ K) im- 
plies a shock velocity of 920 km s ', while the optical filaments 
which reflect shocks traversing denser gas have typical line-of- 
sight velocities of 250 km s"'. Hill et al. (1995) estimated the 
mass of the progenitor of N49 to be about 20 M© based on the 
association with the 10 Myr old association LH 53. A 1720 MHz 
OH maser emission comes from a region toward the southwest. 
N49 is almost certainly the result of a Type II explosion. 

Apparently first detected at mid-IR wavelengths with IRAS 
(Graham et al. 1987), the first mid-IR images of the SNR were 
obtained with Spitzer (Williams, R. et al. 2006), where N49 is 
visible not only in the MIPS 24 and 70 yum bands, but also in the 
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IRAC 3.6, 4.5, 5.6 and 8 //m bands. IRS spectra of the eastern 
limb show strong line emission from O, Ne, Ar, Si, S, and Fe and 
H2 in the 5-37 fj.m band (Williams, R. et al 2006). While lines 
and dust emission compete at the short wavelengths, van Loon 
et al (2010) showed with the low resolution MIPS spectra of the 
bright blob on the eastern limb that dust is the dominant source 
of emission at 70 //m, and estimated a dust mass for N49 of ~0.2 
Me. 

Dust in SNRs can arise either from the SN itself or from the 
ISM. Ejecta dust is thought to play a key role as a coolant in the 
formation of high-redshift galaxies (see, e.g. Morgan & Edwards 
2003). However, given the amount of swept-up material in N49, 
most of the dust in N49 is expected to be of interstellar origin. 

Here we describe observations of N49 obtained with the 
Herschel Space Observatory. Our primary purpose is to explore 
the properties of dust in N49, and to obtain a better estimate 
of the dust mass than was heretofore possible by probing the 
longer wavelength bands, where larger dust grains are expected 
to be more apparent. Combining Herschel data and data obtained 
from other observatories, we construct spectral energy distribu- 
tions (SEDs) of the entire SNR and of a bright blob on the east- 
ern limb, where the SNR is interacting with a molecular cloud. 
(see, e.g. van Loon et al. 2010). We use a simple two-temperature 
model of the SED to estimate dust masses and temperatures. 



2. Observations and data reduction 

N49 was observed with Herschel as part of the science demon- 
stration program associated with the HERITAGE program ( PI: 
M.Meixner). The observations and the data reduction process 
for the PACS green and red bands (100/160 fim) and SPIRE all 
bands (250/350/500 jum) are described by Meixner, et al. (2010). 



gxlO"°4xlO"°8xlO""axlO"° 10"' 



so 



A 

I 



o 

0) 
O 

S3 -50 



g 

3 



IRAC/8 um 



r 

\ E blob ' 



■■■'■■■■'■■ 



Otsuka et al.: Dust in the LMC SNR N49 

2x10"' 4x10"* O.DZ 0.04 0.06 OJMI 5x10"" 0.01 0.01E 0.0Z 

lIIPS/70 um 



lIIPS/24uin 



C 



■■'■■■■'■■■■'■■ 



# 



■■'■■■■'■■■■'■ 




■■''■■ 



i 



SO -50 SO 

0.01 0.02 0.03 0.04 0.05 



-SO 50 

0.15 0.06 



50 SO 

0.15 0.05 



-50 




SO -SO 50 -50 

E <— Relative Position (arcsec) — > W 

Fig. 1. Selected images of the SNR N49. The intensity is in Jy per pixel. The location of the eastern blob (E blob) is indicated by 
the arrow. The radii of the green and yellow circles are 50" and 10" when adopting aperture photometry for the whole SNR and the 
E blob. 



To complement the Herschel data, we used archival 
SpitzerflRSi low- and high-resolution spectra (5-40 yum; PI: 
Gehrz) of the eastern blob (E blob, hereafter), SpitzerllKAC 
3.6/4.5/5.8/8.0 //m and MIPS 24/70 //m images (Meixner et al. 
2006), and the pubUshed 3/7/1 1/15/24 yum fluxes for whole SNR 
(Seok et al. 2008) taken by the AKARIjlRC (Murakami et al. 
2007; Onaka et al. 2007). We indicate the position of the E blob 
in Fig 1 . We rescaled the flux density of the IRS spectrum by a 
factor of 1.25 in order to match the total flux between 20.34 and 
29.46 jim of the IRS spectrum and the MIPS 24 ^im flux. We 
used the reductions of the IRS and the MIPS spectra (52-93 fim) 
of the E blob made by Kemper et al. (2010) as part the SAGE- 
Spec project. The MIPS spectrum has been described in detail 
by van Loon et al. (2010). 

Selected Spitzer and Herschel images are presented in Fig 
1. Although the angular resolution of PACS at 100 and 160 jjm 
is comparable to that of MIPS 24 and 70 yum, the resolution of 
SPIRE is such that it becomes increasingly difficult to distin- 
guish the SNR from the surrounding molecular clouds. The in- 
tensity of the E blob becomes comparable to the southern molec- 
ular clouds at longer wavelengths. There is no obvious differ- 
ence however in the morphology of the SNR as a function of 
wavelength, although at long wavelengths, where the resolution 
is poorer, it becomes increasingly difficult to pick out the shell of 
the SNR. The general interpretation that the SNR is interacting 
with a molecular cloud on the SE is confirmed by the Herschel 
images. 

To estimate flux densities, we performed aperture photome- 
try. We used 50" and 10" radius regions to define the entire SNR 
and the eastern blob, respectively (see Fig 1). We subtracted 
background from an annulus centered on the SNR with inner and 
outer radii of 1 10" and 130". In the IRAC 8 fim image of Fig 1, 
the two regions are indicated as green (50" radius) and yellow 
(10") circles. We performed aperture corrections using observed 
or theoretical point-spread functions. The measured fluxes are 



listed in Table 1 . The absolute flux calibration is 2 and 5 % for 
MIPS 24 and 70 yum bands and is 20 % for PACS (Poglitsch et al. 
2010) and 15 % for SPIRE (Swinyard, Ade, & Baluteau 2010), 
respectively. We adopted these values as the flux density uncer- 
tainty for these bands. We did not include uncertainties associ- 
ated with the subtraction of background, but these are at least as 
large in the SPIRE bands as the uncertainty in the flux calibra- 
tion. For the other bands, the uncertainty in Table 1 corresponds 
to the standard deviation of the background. 

The SEDs are presented in Fig 2 for the entire SNR and for 
the E-blob alone. There are no obvious diff'erences aside from 
normalization in the two spectra. Williams, R. et al. (2006), 
based on an early reduction of the IRS spectrum, argued that 
the emission lines rather than the dust continuum dominate the 
flux in the MIPS 24 /im band. However, our improved reduction 
of the IRS spectrum shows an obvious dust continuum. Our esti- 
mate based on the full 9.7 yum width of the MIPS 24 yum band is 
that emission lines contribute about 38 % of the flux. The fluxes 
plotted at 24 //m in Fig 2 represent the continuum flux, assuming 
that the line contribution is the same for the entire SNR as for the 
region covered by the IRS spectrum. At wavelengths of A ~70 
/im and longer, the radiation from dust grains in the SNR is the 
main luminosity source, van Loon et al. (2010) found dominant 
[O i] 63 fj.m line in the MIPS 70 yum band to contribute only 1 1 % 
of the total emission in that band. Furthermore, based on this, 
the expected [Oi] 146 yum (^Pi-Vq) flux is 3.4-7.2x10"'^ erg 
s"' cm"^ when assuming the electron temperature to be lO'' K 
and a density 10-10^ cm""*. The contribution of [Oi] 146 yum in 
the PACS 160 yum band is -7%. 



3. Dust mass in the N49 system and discussion 

We estimated the dust mass with the line-contribution-subtracted 
data and its temperature through SED fitting. First, we esti- 
mated the contribution of the synchrotron emission to the mid- 
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Table 1. Flux density of N49. (In the last column) (1): AKARI 
(Seok et al. 2008) (2): Spitzer (present work); (3): Herschel 
(present work). 



Table 2. The derived dust mass and temperature. 
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Fig. 2. SED plots of the whole SNR (upper panel) and E blob 
only (lower panel) together with silicate dust fits. The strong 
emission lines of atomic gas and H2 in the IRS spectrum are 
subtracted. The black solid and broken lines are the SEDs of 
warm and cold components, reflectively. The red lines are the 
sum of these two components. 



to far-IR region. Dickel & Milne (1998) measured the flux den- 
sities at 3.5/6/13/20 cm for the whole SNR. By extrapolating 
the flux densities at this radio wavelength to the mid-IR region, 
we found that the contribution is negligibly small (<2 %) to the 
24/70/100/160/250 jum flux densities. In the SPIRE 350/500 yum 
bands, the contribution is large; ~14 % for 350 yum and -79 % 
for 500 /^m. Before fitting, we subtracted the contribution of the 
synchrotron emission to the 350 yum flux density for the whole 
SNR and E blob, assuming that the contribution can be also ap- 
plied to the E blob. In the fittings we excluded 500 fiva data. 
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When we assume optically thin thermal radiation at a dust 
temperature T^j of a component / and a single dust species with 
the radius a, the observed flux density Fy is written by 



F, = {^apoA X ^ MjjB,(Tdj)Qyj, 



(1) 



where M,/ ,■ is the dust mass of the component /, By(Tdj) is the 
Planck function of Tdj, and Qyj is the absorption efficiency, p is 
the dust grain density and D is a distance to the observed dust 
from us (Kwok 2007). Here, we adopted D=50 kpc (Schaefer 
2008). We considered two types of grains, amorphous carbon 
and silicates. For amorphous carbon the optical constants were 
taken from Zubko et al. (1996). We adopted p=2.26 g cm"^. To 
simplify the model, we assume a = 0. 1/vm uniform grain radius. 
Next, we considered astronomical silicate only. The optical con- 
stants were taken from Draine & Lee (1984) and a p of 3.6 g 
cm-^ was adopted. Due to the lower absorption efficiency at >70 
pm than amorphous carbon, a larger dust mass would be esti- 
mated. The dust mass estimated for the silicate-only case would 
be an upper limit. In the fittings we found the best fit M,/,, and 
Tdj to the observations by least-square technique. In this pro- 
cess, Mrf , and Tdj are the free parameters, and all the data are 
equal of weight. 

We first considered a single temperature model to fit the 
SED. For the whole SNR, the best fit had a T^ of 33 K and 34 K 
for amorphous carbon and for silicate respectively, and a derived 
Md of 3.6 Mo and 4.6 Mq. For the E blob, Td is 35 K (carbon 
case) and 39 K (silicate case), and Md is 0.47 Mq and 0.60 Mq. 
The temperature is fairly similar to the temperature derived by 
Graham et al. (1987) from IRAS data. However, we found that a 
better fit could be obtained if we considered a model with two- 
temperatures instead of one. The one-component fittings cannot 
simultaneously explain the <24 pm and >160 pm flux. The im- 
plied dust masses and temperatures for the two-temperate model 
are summarized in Table 2 for the two types of grains. 

The best-fit two-temperature model, assuming silicate 
grains, for the SEDs is shown in Fig 2. The warm component 
has a temperature of ~60 K, and is surely due to dust in shocked 
gas, while the cold component has a temperature of ~20 K, sim- 
ilar to that found in molecular clouds (16-24 K, Bernard et al. 
2008; Gordon et al. in 2010). The contribution from the cold 
component to the total flux in the Herschel bands is >80 %. 
The dust mass and temperature of the warm component in the 
E blob are very close to van Loon et al. (2010; 0.2 M© and 43 
K). The total mass of dust associated with N49 is about 10 M© 
and is dominated by dust at cold temperatures, dust which was 
difficult to observe with Spitzer. It is tempting therefore to iden- 
tify the warm component with dust in the X-ray plasma and the 
cool component not with the SNR, but with the molecular cloud. 
However, this assignment is very uncertain, because dust tem- 
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Fig. 3. The relation between the dust mass and age. The data are 
taken from Seok et al. (2008) for 0509-67.5, 0519-69.0, N132D, 
and N49B; WilUams, B. et al. (2008); Sandstrom et al. (2009) 
and Rho et al. (2009) for IE 0102.2-7219; WilUams, B. et al. 
(2006) for N132D, N23, 0453-68.5, and N49B; van Loon et al. 
(2010) and present work for N49. 

perature variations are also expected as a result of the large vari- 
ations of post-shock gas density that are known to exist in N49. 

To place SNR N49's dust mass estimate into the context of 
SNR development, we plot the dust mass vs. age for all LMC 
and SMC SNRs studied to data, as noted in Fig 3. Note that 
except for N49 the dust masses are estimated based on <70 ycca. 
data. The relation can be with the free expansion phase early on, 
followed by a Sedov phase (radius oc (age)^^^; cf. Badenes et al. 
2010). 

The mass of the warm dust component in N49 is 0. 1 Mq 
assuming amorphous carbon and 0.4 M© assuming silicates. 
Adopting a dust-to-gas mass ratio of ~250 (Meixner et al. 2010) 
and assuming that only the warm dust is associated with shocked 
gas, the swept-up mass of the SNR ranges from 25 to 100 Mq. 
This is comparable to the mass estimated from X-ray observa- 
tions, 200 Mq (Hughes et al. 1998), assuming the true grain 
characteristics are similar to those of silicates. A much higher 
(25 X larger) swept-up mass is estimated if we included the cold 
dust in our estimate, which suggests that in the context of our 
simple model there is a significant molecular cloud component 
to the cold dust we associate with the SNR, and that the shocks 
propagating into the molecular cloud may be too slow to pro- 
duce X-rays. While more sophisticated modeling of N49 may 
refine the mass estimates somewhat, the basic limitation in ar- 
riving at a more precise picture of the nature of the interaction 
with the molecular cloud and the amount of shocked molecular 
gas in N49 is the limited angular resolution of He rschel at longer 
wavelengths. 
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